SHEP 02-17 
hep-ph/0207260 

Gamma ray bursts as probes of neutrino mass, quantum 

gravity and dark energy 



Sandhya Choubey and S.F. King 

Department of Physics and Astronomy, University of Southampton, 
Highfield, Southampton SOI 7 IB, J, UK 



Abstract 

We calculate the time delays of neutrinos emitted in gamma ray bursts due to the effects 
of neutrino mass and quantum gravity using a time dependent Hubble constant which can 
significantly change the naive results presented hitherto in the literature for large redshifts, 
and gives some sensitivity to the details of dark energy. We show that the effects of neutrino 
mass, quantum gravity and dark energy may be disentangled by using low energy neutrinos 
to study neutrino mass, high energy neutrinos to study quantum gravity, and large redshifts 
to study dark energy. From low energy neutrinos one may obtain direct limits on neutrino 
masses of order 10 -3 eV, and distinguish a neutrino mass hierarchy from an inverted mass 
hierarchy. From ultra-high energy neutrinos the sensitivity to the scale of quantum gravity 
can be pushed up to Eqq ~ 5 x 10 30 GeV. By studying neutrinos from GRBs at large redshifts 
a cosmological constant could be distinguished from quintessence. 
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1 Introduction 



Gamma ray bursts (GRBs) are amongst the most distant, energetic and enigmatic astrophysical 
phenomena known. Understanding GRBs is arguably the most outstanding question in astronomy, 
and one which may be answered by a plethora of gamma ray observatories such as INTEGRAL, 
SWIFT and BATSE, and corresponding infrared and optical telescopes such as REM and LT [1] . 
It is well known that 99% of the energy of a supernova is emitted in the form of neutrinos, and 
therefore it is widely expected that GRBs are similarly a copious source of neutrinos which may 
be detected in future neutrino telescopes [2, 3, 4, 5]. Within this decade it is therefore likely that 
GRBs will become much better understood, and their exact nature and mechanisms which drive 
their internal engines will be revealed. For example it may turn out that a GRB results from 
the core collapse of a very massive supernova to a compact rotating black hole with the energy 
emitted in beamed relativistic fireball jets containing copious neutrino fluxes [6, 7, 8, 9, 10, 11] 
Alternatively the GRB engine could result in the emission of beamed earth-sized cannonballs [12]. 

In this paper we are not concerned with detailed models of GRBs, but instead regard them as 
a high intensity, high energy neutrino beam with a cosmological baseline. We shall be interested 
in the time delay of the arrival of neutrinos. We show that the time delay may be used as a 
probe of three physical effects: (i) neutrino mass, (ii) quantum gravity, (iii) dark energy. The 
time delay due to neutrino mass has been noted earlier in [2] while effects of quantum gravity 
on the time of flight have been considered in [13, 14, 15, 16] for high energy photons and in 
[17, 18] for neutrinos. However none of the above papers considers the time delays due to massive 
neutrinos in the presence of quantum gravity, although [17] gives the dispersion relation for this 
case. Moreover no studies to date have calculated time delays using a formalism which correctly 
takes into account the time dependence of the Hubble constant due to matter and dark energy. In 
our study we consider the time dependence of the Hubble constant due to matter and dark energy 
and show that they can change the naive results by more than 100% for z > 1. For large redshifts 
the results for time delays due to neutrino mass and quantum gravity are sensitive to the nature 
of dark energy. We show that the three effects may be disentangled by using low energy neutrinos 
to study neutrino mass, high energy neutrinos to study quantum gravity, and using large redshifts 
to study dark energy, leading to the results stated in the abstract. 

For the determination of the neutrino mass, in principle one could compare the arrival time 
of the massive neutrinos with the arrival times of the photons emitted in the GRB, assuming 
them to be emitted at the same time. However this would be plausible only if the GRB were a 
point source in vacuo, which it is not. In any realistic GRB model the photons are trapped inside 
the fireball and are released much later - the exact amount of time delay being highly model 
dependent. Another strategy might be to compare the arrival times of the low energy neutrinos 
with that of the ultra-high energy ones, since (as we will show) the ultra-high energy neutrinos 
suffer negligible time delay due to mass. But again one would expect the low energy neutrinos to 
be produced thermally, as in supernovae, leading to a model dependent time delay. 

However there are alternative strategies which could overcome these problems. To begin with, 
if the neutrinos are hierarchical in mass, m\ <C m 2 <C m 3 , then due to their mixing, neutrinos 
of the same energy will arrive at the detector in three bunches, corresponding to the three mass 
eigenstates, the arrival time only depending on their mass. We can then compare the arrival times 
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of the different neutrino mass eigenstates and put limits on the neutrino mass. This represents 
a "clean strategy" and indeed similar arguments for constraining neutrino mass from time de- 
lays using supernova neutrinos have been used in the literature (see [19] and references therein). 
However we will show that with GRB neutrinos the sensitivities possible are better by orders of 
magnitude. 

For constraining models of quantum gravity and dark energy we can use arrival times of the 
ultra-high energy neutrinos, which are produced over many decades of neutrino energy. We can 
compare the arrival times of the high energy neutrinos of different energies - which unlike the low 
energy thermal neutrinos are expected to be produced at almost the same time at the source - and 
put limits from such observations. Thus, the results in this paper do not rely on the comparison of 
the arrival times between photons and neutrinos, which would involve the uncertainties discussed 
above. However for simplicity, we shall calculate arrival time delays relative to a hypothetical low 
energy photon which is assumed to be emitted at the same time. This is for convenience only; in 
a realistic search strategy what will be important will be the comparison of time delays between 
neutrinos of either the same energy, or between high energy neutrinos of two different energies as 
discussed above. 

2 Formalism 

For a neutrino of mass m and energy E(t), the momentum p(t) in the presence of the effects of 
quantum gravity (QG) with an effective energy scale Eqc is given to lowest order as [17] 



where £ = ±1. It is worth briefly discussing the origin and reliability of the QG corrections to the 
dispersion relation in Eq.l. Ellis et al [16, 18] have shown that quantum fluctuations in space-time 
foam lead to a energy dependent perturbed gravitational background metric 



where \u/c\ ~ E/Eqc, Eqc being the scale at which the quantum gravity effects set in. Such a 
change in the metric violates Lorentz invariance at some high scale Eqc and changes the dispersion 
relations of the particles. The linear energy dependence of such Lorentz invariance violating (LIV) 
terms in these models come mainly from the gravitational recoil effects. An identical form for 
the dispersion relation Eq.l is obtained by studying massive spin-1/2 fields in the framework of 
loop quantum gravity where the LIV breaking correction term in Eq.l arises due to the discrete 
structure of space-time at Planck scale [17]. 

In fact the last term in Eq.l can arise in any theory which violates Lorentz invariance. The 
energy dependent form of the corrections may arise from a dimension 5 LIV operator of the 
form c^xip^^D^Dxip. The dimension 4 operator c^ip^^D^tp would lead to an energy independent 
correction to the dispersion relation while the dimension 6 operator cipj^D^D 2 ^ would lead to 
only tiny energy independent corrections. Hence one may regard the dispersion relation in Eq.l 




(1) 



9ij — 900 



1, 9oi = Ui/c 



(2) 
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as a consequence of a generic type of LIV which may arise in some particular QG theories. From 
the point of view of measuring neutrino masses, the QG corrections to the dispersion relation 
represent a possible effect which may threaten to swamp the neutrino mass effect, and therefore 
it is important to include the largest imaginable such effects as we do here. 

Assuming Eq.l the time dependent neutrino velocity v(t) is then given to lowest order in terms 
of the observed neutrino energy E as 1 

dE(t) „ ( el a\t) 3 a(t )\ 

where we have expanded in Cq <C 1 and e[, C 1 where 

e = -^, e = — ■ (4) 

&QG 

The time-dependence of the neutrino velocity v (t) arises from the expansion of the universe which 
redshifts the neutrino de Broglie wavelength A(t), and reduces their momentum and hence their 
velocity. The neutrino momentum p(t) is related to the cosmic scale factor a(t) through 

P(t) = Ago) = ago) 

p(t ) X(t) a(t) { ' 

Now suppose that low energy photons 2 and high energy neutrinos are emitted from a GRB 
source at time t e and the low energy photons arrive on Earth at time t while the neutrinos arrive 
at time t v . The low energy photons will travel at the speed of light c, while the high energy 
neutrinos will travel with time dependent velocity v (t) which may be smaller than c due to the 
finite neutrino mass, and may be smaller or larger than c due to the dispersive effects of quantum 
gravity. The co-moving distance x between the GRB source and the Earth calculated in terms of 
the low energy photons is given by 3 



x = / ~77\ W 
Jt e a(t) 



where a{t) is the cosmic scale factor. The co-moving distance x between the GRB source and the 
Earth calculated using the neutrinos which travel freely with time dependent velocity v (t) is 



X 



v(t)dt 
a(t) 

/■*> v{t)dt + v(t)dt 
Jt e a(t) Jto a(t) 

1 Strictly E is the neutrino energy at the time that the photons are observed on Earth E = E(to), but since 
the neutrinos arrive a short time later this is to excellent approximation equal to the observed neutrino energy. 

2 As discussed in the Introduction, the time delay is defined in terms of low energy photons for convenience only. 
What will matter in a practical search strategy is the comparison of different time delays between neutrinos of 
cither the same energy, or between high energy neutrinos of two different energies. 

3 While we have used he low energy photons for defining the co-moving distance, Eq.6 holds even for low energy 
massless neutrinos which travel at speed c. 
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The co-moving distance in Eq.7 must be equal to that calculated using the low energy photons or 
low energy massless neutrinos in Eq.6. 

If we now assume that the neutrino velocity and the cosmic scale factor do not change much 
over the time scale i — 1„, then equating Eq.7 to Eq.6 and using Eq.3 we find the time delay, to 
leading order in and e , 



At 



2 Jt e alin) 2 Jt e a z (t) 



a(i ) 2^ u A e a 2 (i) 



The time delay in Eq.8 may be expressed in terms of an integral over redshift defined as 1 + z' = 
a(t )/a(t), and Hubble constant defined as H = a/a, 

Ai « |/ 2 - l&I-u (9) 
where I n are "redshift moments" of the inverse Hubble constant, 

In = So (l + z')"H(z'Y (10) 

It is interesting to compare the expression for the time delay in Eq.9 to the result for the total 
time taken T for the low energy photon or massless neutrino to travel from the GRB source to 
the Earth, 4 

T = t -t e = I u (11) 

It is also interesting to compare to the result for the proper distance D of the GRB source to the 
Earth at the photon arrival time, 5 

D = a(t )x = cl . (12) 

The upper limit of the integrals in Eq.10 z represents the redshift of the GRB source, where 
1 + z — a(t )/a(t e ), and the Hubble constant H(z) for a flat universe is given by 

H 2 (z)/H^ = n M (l + zf + n DE e I(z) (13) 

where 

I(z) = 3 f\l + w(z'))-^- (14) 
Jo 1 + z 

where the equation of state for the dark energy w, defined as the ratio of its pressure to its density 
w = p/p, in quintessence models has a redshift dependence. For constant w we have 

H\z)/Hl = Q M (1 + zf + n DE (l + zf^ (15) 



4 The time T represents how far back in time the GRB took place. As z — > oo, T — > To where To is the age of 
the universe. 

5 D represents the actual distance of the GRB from the Earth as measured now. As z — > oo, D — > Hu where 
Hd is the horizon distance that light could have travelled since the beginning of the universe, which represents the 
size of the observable universe. 
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Figure 1: The integrals /„ in units of the Hubble time t H rs 14 billion years as a function of the 
red shift z for three different cosmological models: CONST (solid), SUGRA quintessence (dots), 
INV quintessence (dashes). Iq is the distance of the GRB in light years, ii is the date (years B.C.) 
that the GRB exploded, I2 is used to calculate the time delay At due to neutrino mass, and i_i 
is used to calculate the time difference due to quantum gravity effects as discussed in the text. 

and for a cosmological constant A with w — — 1 and Vl DE = f2 A this reduces to, 

H 2 (z)/H 2 =Q M (l + z) 3 + Q A . (16) 

The present day ratio of matter (M) density to critical density is Qm ~ 0.3, the present day 
ratio of dark energy (DE) density to critical density is VLde ~ 0.7 and the Hubble constant is 
H « TOkms^Mpc" 1 . 

We shall consider three different cosmological models which were recently parametrised in 
[20]: cosmological constant (CONST) with potential V ~ A 4 ; supergravity (SUGRA) inspired 
quintessence field (Q) model with potential V(Q) ~ l/Q u e Q / 2 ; inverse (INV) power law quintessence 
field (Q) model with potential V(Q) ~ 1/Q 6 - The equation of state w(z) for the SUGRA and 
INV models does not have a simple analytic form, but it may be parametrised as discussed in [20] , 
and we shall use the parametrisation given there in this paper. Each of these models has w(z) 
which varies with z in such a way as to lead to a "tracking" behavior. This implies in particular 
that the value of the dark energy density remains of the same order as the matter density for 
large redshifts. By comparison the ratio of the cosmological constant energy density to the matter 
density falls as 1/(1 + z) 3 and rapidly becomes negligible at large redshifts. 

In Figure 1 we calculate /„ in Eq.10 for n = —1, 0, 1, 2 in units of the Hubble time tn — V-Hio ~ 
14 Gyr as a function of redshift z for the three different cosmological models CONST, SUGRA, 
INV defined above. The physically relevant quantities D, T and At are simply related to J , h, 
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I 2 and 7_i as shown in Eq.12,11,9. Figure 1 shows that as n increases /„ becomes increasingly 
insensitive to redshift (flatter curves) and to the particular cosmological model (closer spacing 
between solid, dotted and dashed curves) at large redshift. Both effects can be simply understood 
as being due to the factor of (l + z') n in the denominator which tends to suppress the contributions 
from the higher redshift parts of the integration region for larger n. The difference between the 
curves for I and those for I 2 in Figure 1 represents the error that would be made if the time 
delay were calculated using the naive formula 



Af 





(17) 



which ignores the effect of neutrino redshift as has been done up till now in the literature rather 
than the correct formula for At in Eq.9. For z = 1 the error incurred by using the naive formula 
can clearly be seen to be of order 100%, with the error rapidly growing for larger redshifts. It 
is also clear from Figure 1 that there is some sensitivity to the nature of dark energy for large 
redshifts. 

Finally it is interesting to consider the matter dominated limit of our time delay result corre- 
sponding to Qm = 1, &de = 0. In this limit the expression for the time delay in Eq.9 reduces 
to 



At 



matter 



5#n 



1 - 



1 



(1 + zfl 2 



Ha 



(1 + *) 



1/2 _ 



(18) 



Time delays due to quantum gravity effects in the matter dominated limit were also considered 
in [16]. However the result quoted there corresponds to the second term in Eq.17 rather than the 
second term in Eq.18 which correctly takes into account the expansion of the universe. 



3 Results 

The results in this section are based on the full formula for At in Eq.9, using the Hubble constant 
for a flat universe calculated using Eq.13 for the different dark energy models. 

Figure 2 shows the time delays of the neutrinos against observed neutrino energy E for a 
fixed GRB redshift of z — 1 and assuming the cosmological constant model. In each panel the 
downward sloping dotted lines give the time delay due to the effect of neutrino mass which in 
the upper panels is chosen to be m = 0.05 eV corresponding to the "atmospheric neutrino mass" 
defined as the square root of the atmospheric mass squared splitting [21], and in the lower panels 
we take m = 0.005 eV corresponding to the "solar neutrino mass" which is the square root of the 
large mixing angle solar mass squared splitting [22, 23]. To probe neutrino mass corresponding to 
the solar scale the detector has to observe neutrinos with E ~ few 10 MeV with time sensitivity 
of a few milliseconds. The time delays corresponding to the atmospheric scale are higher and 
should be easier to detect in the planned neutrino telescopes. 

The upward sloping dashed lines show the gravitationally induced time delay of the neutrinos. 
The upper panels give the time delay when the quantum gravity energy scale corresponds to the 
Planck scale. The lower panels are for the case where the quantum gravity corrections become 
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important at 10 22 GeV. The left panels have £ = +1 corresponding to a negative time delay due 
to quantum gravity which tends to cancel with the positive time delay from the effect of neutrino 
mass. For this case the higher energy neutrinos arrive earlier. In the figure we have plotted 
the absolute value of the difference between their arrival times. The right panels have £ = — 1 
corresponding to a positive time delay due to quantum gravity which reinforces the positive time 
delay from the effect of neutrino mass. For both the cases we observe that the quantum gavity 
effects become very important for neutrinos arriving with energies greater than a few GeV. The 
ultra high energy neutrinos travelling cosmo logical distances can put severe constraints on Eq G . 
Note that the arguments in [18] suggest that fermions travel more slowly than c and hence that 



The solid lines in the figure show the time delay of the neutrinos due to the sum of the dotted 
and dashed lines when both the effects are present. For the case where £ = +1 the two effects can 
cancel each other when from Eq.9 we have At = which occurs at an energy 



For £ = — 1 even though the two effects reinforce each other, we see from Eq.9 and fig. 2 that the 
minima in At comes at exactly the same energy. 

Figure 2 shows that time delay due to neutrino mass is important for the lower energy neutrinos 
with Eq < E' Q , while the time delay due to quantum gravity effects is important for the higher 
energy neutrinos with E > E' . This enables the two effects to be disentangled and treated 
separately. 

3.1 Neutrino Mass Limits from Low Energy Neutrinos 

In this section we shall focus on neutrino events with energy E < E' in order to put limits on m. 
This implies that neutrinos with energies less than ~ 100 MeV can be effectively used to study 
their mass. Thermal neutrinos with energies between 10 — 100 MeV are expected to be produced 
in GRBs [2, 11]. Thermal neutrinos, coming from GRBs which are energetic enough, which are 
not very far away and which are probably beamed towards the Earth, should be detectable in the 
future km 2 ice detectors like IceCube. These lowest energy neutrinos are expected to be detected 
mainly via P e + p — > n + e + where the positrons lead to an increase in the low PMT noise. This 
detection method which forms part of a Supernova Early Warning System may also be used to 
detect low energy neutrinos from GRB's [2, 11]. However the determination of neutrino mass 
using the time delay techniques requires sensitivity to both arrival times as well as energy of the 
incoming neutrino. For a detector like IceCube since the inter-PMT distance is large the energy 
resolution is expected to be poor, particularly for the lower energy neutrino events. However 
energy resolution of the proposed megaton Super-Kamiokande type water Cerenkov detectors 
such as Hyper-K, UNO and TITAND should be good and can be effectively used for putting 
direct limits on neutrino mass. We next calculate the number of events that these detectors would 
observe from a GRB event. 

One can estimate the number of thermal neutrinos emitted in a typical GRB by the following 
argument [2, 11]. Assume that the energy of the GRB in photons is E 1 and the radius of the GRB 




(19) 
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Figure 2: Time delay At due to the effect of both neutrino mass and quantum gravity against 
observed neutrino energy E for z ~ 1 for the cosmological constant model. 
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Figure 3: Number of events expected from thermal neutrinos in a one megaton water cerenkov 
detector as a function of the GRB redshift z. 
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fireball is R, then the photon temperature T 7 is given by T 7 = (2E 1 /aVh 1 ) 1 ^ 4: , where V = ^irR 3 , 
a is the Stefan-Boltzmann constant and /i 7 = 2 is the number of degrees of freedom. Similarly 
the neutrino temperature T u is given by T v = (2E u /aVh l/ ) 1 ^ 4 , where the number of degrees of 
freedom for the neutrinos is h u = 2 x 3 x | and E u is the energy emitted in neutrinos. We shall 
assume E v ~ 100-E 7 [2, 11]. The total number of neutrinos emitted is given by N v = E V /(E V ), 
where (E u ) = 3.1571, is the average energy of the neutrinos, assuming a Fermi-Dirac distribution. 
From this we obtain 

/ P \ 3/4 / R \ 3 / 4 
Ar, = 7 xl0 57 — ^— (-^—) ■ (20) 
\10 54 erg S y V 100km/ v ; 

The neutrino flux at Earth is 

$ = (21) 

V 47T J D 2 fi beam 1 ; 

where D is the distance of the GRB given by Eq.12 and fibeam is the beaming solid angle of the 
GRB corresponding to a cone of opening angle #beam, which we assume to be of order 6*beam ~ 1° 
[24]. The total number of events observed at the detector by v e capture on protons is given by 

N = N P J <f)(E)a(E)dE (22) 

where the v e flux <f>(E) = ^> v f FD (E), f FD (E) being the Fermi-Dirac distribution function and 
we assume that the total neutrino flux <& u is evenly distributed in all the six species. The reaction 
cross-section <r(E) for u e + p — > n + e + is given by 



a(E) = 0.94 x 10" 43 (E - Q) [{E - Qf - 



ml) V2 cm 2 (23) 



where Q = 1.293 MeV is the neutron-proton mass difference and m e is the electron mass. N p = 
j^NaMd is the number of free protons in water, where Na is the Avogadro's number and Mp is 
the fiducial mass of the detector. The integration in Eq.22 is done over observed neutrino energies. 
In Figure 3 we show the number of events expected in a one megaton water Cerenkov detector 
from thermal neutrinos produced in a GRB, which is beamed towards us with a beaming angle of 
1° [24]. The number of expected events are seen to be large for the GRBs with low redshift. 

In Figure 4 we show the time delay At in seconds as a function of eo = mc 2 / E$ for redshift 
z = 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 10, for the cosmological constant model. In this figure we assume 
that the neutrino energy satisfies E < E' so that the effects of quantum gravity can be neglected 
and hence we can set = 0. For a fixed redshift the time delay increases as el = m 2 c A /El giving 
larger time delays for more massive neutrinos and lower observed energies, as seen in Figure 2. We 
are interested in measuring the smallest neutrino masses, which for a fixed eo corresponds to the 
smallest observed energies E and the smallest time delays. Since the neutrinos are released from 
the GRB in bursts on the time scale of milliseconds, the smallest time delays that are meaningful 
will also be milliseconds and from the figure this corresponds to e = mc 2 /E ~ 10" 10 . For the 
lowest detectable energy neutrinos of E ~ lOMeV, e ~ 10~ 10 corresponds to neutrino masses of 
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Figure 4: The time delay At in seconds as a function of eo = mc 2 / Eq for redshift z = 
0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 10, assuming the cosmological constant model and taking e' = 0. 




Figure 5: Variation in time delay At in the millisecond region with eo = mc 2 / Eq for red shift of 
z = 0.01,0.02,0.05,0.1,0.2,0.5, 1, 10, for the three different cosmological models (taking e = 0): 
CONST (solid), SUGRA quintessence (dots), INV quintessence (dashes). 
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m ~ 10~ 3 eV. In order to discuss the limits on neutrino masses in more detail we must focus on 
the millisecond region of this plot. 

In Figure 5 we show a blow-up of Figure 4 in the important millisecond region. We have also 
included the effect of different cosmological models in this figure, which can be important for high 
redshifts. From Figure 5 we see that time delays of 10~ 3 s correspond to eo = 7 x 10 -10 — 8 x 10~ n , 
over a range of redshift z = 0.01 — 10. For the lowest energy neutrinos of Eq ~ lOMeV this range 
of e corresponds to a range of neutrino masses of m = 7 x 10~ 3 — 8 x 10~ 4 eV. Although the 
smallest neutrino masses are associated with the highest redshifts the number of events expected 
in the detector falls sharply with increasing redshift. Since i_i is not very sensitive to z beyond 
z ~ 1, the fall in the number of events with distance is more acute than the rise in the time 
delay, and hence GRBs at lower values of redshift would be better suited for determination of the 
neutrino mass. 

In addition to being sensitive to the absolute masses of all the three neutrino states, GRB 
neutrinos have the potential to probe the mass hierarchy as well. Though from solar neutrino 
data we know that the sign of Am^ = Am,Q is positive [22, 23] (Am| = mf — m 2 ), there is still 
an ambiguity in the sign of Am| 2 = Am„ (m . It would be hard to determine the sign of Am^ in any 
of the current and planned long baseline oscillation experiments and only a neutrino factory would 
be able to resolve this ambiguity. The delay in arrival times for the neutrinos can in principle be 
used to determine the neutrino mass hierarchy. Neutrinos arrive in three bunches corresponding 
to the three mass eigenstates. The neutrinos are detected via the electron antineutrino reaction, 
so for each mass eigenstate we must multiply the detection rate by the probability \U ei \ 2 that the 
mass eigenstate corresponding to the mass m 8 contains an electron antineutrino. Since \U e3 \ 2 is 
small m 3 will have the smallest component of electron antineutrino and hence smallest number of 
events at the detector. In practice the states with mass mi and vri2 may be most easily detectable 
since \U e \\ 2 and \U e 2\ 2 are not too small. The heaviest mass eigenstate will arrive last and have the 
largest time delay. In hierarchical models the heaviest mass is m 3 while in models with inverted 
mass hierarchy m 3 is the lightest. So depending on whether fewer events are detected for the 
neutrinos which arrive last (earliest) it would in principle be possible to conclude that the mass 
hierarchy is normal (inverted). 

3.2 Quantum Gravity Limits from High Energy Neutrinos 

The dispersion of velocity due to quantum gravity effects leading to delay in arrival times of the 
neutrinos is most significant at the higher energy end of the GRB neutrino spectrum, as evident 
from Figure 2. In fact for a given neutrino mass m, all neutrinos with energy Eq > E' Q can be used 
to put limits on the scale of quantum gravity Eq G . Since neutrinos with a wide range of energies 
from 1 — 10 9 GeV are expected to be produced in GRB fireball [6, 7, 8, 9, 10], the detection of their 
time delays is the most powerful method of studying/constraining models which predict dispersion 
relations given by Eq.l. 

In this section we therefore focus on neutrinos with energy Eq > E' in order to set limits on 
the quantum gravity scale. Figure 6 shows the time delay expected as a function of e' = E /Eqq 
for redshifts of 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1 and 10, assuming the cosmological constant model 
and also that E > E' so that we can set e = 0. As also seen in Figure 2 the time delays in Figure 
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Figure 6: The time delay At in seconds as a function of e' = Eq/Eqg for redshift z = 
0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 10, assuming the cosmological constant model and taking e = 0. 




Figure 7: Variation in the quantum gravity induced time delay At in the millisecond region with 
e' = E /Eq G for red shift of z = 0.01,0.02,0.05,0.1,0.2,0.5, 1, 10, for the three different cosmo- 
logical models (taking e = 0): CONST (solid), SUGRA quintessence (dots), INV quintessence 
(dashes). 
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6 increase linearly with energy and hence with e' . The time delay also increases significantly with 
the redshift z as i_i has a sharp z dependence. Since we want to restrict ourselves to minimum 
time delays of a millisecond, which corresponds to the variability time of the GRBs, we show in 
Figure 7 the blow up of the region with At ~ 10~ 3 seconds. We also show the effect of the different 
cosmological models on time delays on the same plot. For z = 0.01 — 10 the range of e' which 
would give At < 10~ 3 seconds is e' < 1.5 x 10~ 19 — 2 x 10 -22 . For the highest energy neutrinos 
with E ~ 10 9 GeV expected from the GRB, this would translate into E QG > 6 x 10 28 - 5 x 10 30 
GeV. This should be compared with the bounds set by the time delays of the photons observed 
from GRBs [13] where the tightest limit obtained from the highest energy observed events is 
Eqg > 8.3 x 10 16 GeV [15]. In the extreme ultra-relativistic regime the effect of neutrino mass is 
negligible and the dispersion of the neutrinos in the quantum space-time foam would be almost 
the same as that of photons, the latter having the advantage of being easier to detect. However 
ultra-high energy photons are subject to the GZK cut-off and beyond that for better bounds on 
Eqg one can use ultra-high energy neutrinos. Bounds on Eqg from non-observation of dispersion 
effects of the quantum space-time foam in neutrino oscillations experiments is slightly stronger 
with E QG £ 10 22 GeV [25]. This upper limit may be improved in the forthcoming long baseline 
experiments using neutrino superbeams. But it should be noted that this effect will show up in 
neutrino oscillations only if the dispersion due to quantum gravity is flavor dependent. However 
the method using time delays will be sensitive even if the effect is flavor independent. Thus the 
observation of time delays of GRB neutrinos is the most promising way of phenomenologically 
testing quantum gravity. 

3.3 Dark Energy Limits from High Redshift Events 

Different dark energy models give different values for the Hubble constant H(z) at a given redshift. 
This can have significant impact on the time delays. Since the difference in H(z) is maximum at 
the highest redshifts, we have to look at time delays of neutrinos coming from the farthest GRB 
events to probe dark energy. Since from Figure 1 we see that I-i has a much better sensitivity 
to dark energy than I2, time delays arising from quantum gravity would be a better probe for 
the different cosmological models. Hence it is the ultra high energy neutrinos which would be 
potentially more sensitive. In Figure 8 we show a further blow-up of the time delay plot due to 
quantum gravity effects, with At between 5 x 10~ 3 — 10~ 2 seconds and e' between 10~ 21 — 10~ 20 . 
We show the plots for z =1, 2, 4, 6 and 10. The three different line types correspond to the 
three different models for the dark energy that we have considered in this paper. The time delay 
obtained depends on the redshift as well as the cosmological model. For z > 3 the time delays 
for different redshifts and different cosmologies begin to overlap. This ambiguity can be resolved 
by using smaller redshift events to determine Eqg, since for such events the effect of dark energy 
is reduced. Once Eq G is determined then one can use the high redshift GRB neutrinos to probe 
dark energy. In order to successfully probe the correct cosmological model it is important that the 
redshift of the distant GRBs can be determined accurately enough from their afterglow and that 
the scale of Eqg can be ascertained well enough by looking at the time delay for neutrinos from 
GRB from lower redshifts. If this can be achieved then information on At for neutrinos coming 
from the very distant sources can be used with information on z and the scale of Eq G to probe 
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Figure 8: Variation in the quantum gravity induced time delay At in the millisecond region with 
e' = E /Eq G for high red shifts of z — 1, 2, 4, 6, 10 (from right to left) for the three different cos- 
mological models (taking e = 0): CONST (solid), SUGRA quintessence (dots), INV quintessence 
(dashes). 



the correct cosmological model. 



4 Discussion and Conclusion 

In this paper we have calculated the time delays of neutrinos emitted in gamma ray bursts due to 
the effects of neutrino mass and quantum gravity. Our results are based on the formula for At in 
Eq.9, using the Hubble constant for a flat universe calculated using Eq.13 for different dark energy 
models. This formalism correctly takes into account the time dependence of the Hubble constant 
due to matter and dark energy, and can change the naive results in the literature by more than 
100% for z > 1. We have shown that the effects of neutrino mass, quantum gravity and dark energy 
may be disentangled by using low energy neutrinos to study neutrino mass, high energy neutrinos 
to study quantum gravity, and large redshifts to study dark energy. From low energy neutrinos 
one may obtain direct limits on neutrino masses of order 10 -3 eV, and distinguish a neutrino mass 
hierarchy from an inverted mass hierarchy. From ultra-high energy neutrinos the sensitivity to the 
scale of quantum gravity can be pushed up to Eq G ~ 5 x 10 30 GeV. By studying neutrinos from 
GRBs at large redshifts a cosmological constant could be distinguished from quintessence. 

For convenience we have calculated all time delays with respect to a hypothetical low energy 
photon, assumed to be emitted at the same time from a point source as the neutrino of a given 
energy. We emphasise that what is important in practical search strategies is not these time delays 
themselves, which will be unmeasurable due to the uncertainties in the emission characteristics 
of low energy photons, but rather the comparison of time delays between neutrinos of either the 
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same energy, or between high energy neutrinos of two different energies. As already mentioned, 
if the neutrinos are hierarchical in mass then neutrinos of the same energy will arrive at the 
detector in three bunches, and we can then compare the arrival times of the different neutrino 
mass eigenstates and put limits on the neutrino mass, as in the case of supernova neutrinos [19] 
but with sensitivities better by orders of magnitude. For constraining models of quantum gravity 
and dark energy one can use arrival times of the ultra-high energy neutrinos, and can compare 
the arrival times of the high energy neutrinos of different energies. 
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